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Detection of HD in emission towards Sagittarius B2 
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Abstract. The giant molecular cloud Sagittarius B2, located near the Galactic Centre, has been observed in the 
far-infrared by the ISO Long Wavelength Spectrometer. Wavelengths in the range 47-196 fim were covered with 
the high resolution Fabry-Perot spectrometer, giving a spectral resolution of 30-40 km s~^. The J = 1 ^ and 
J = 2 ^ 1 rotational transitions of HD fall within this range at 112 /j,m and 56 ^m. A probable detection was 
made of the ground state J = 1 ^ line in emission but the J = 2 — > 1 line was not detected above the noise. 
This allowed us to calculate an upper limit on the temperature in the emitting region of approximately 80 K and 
a value for the deuterium abundance in the Sgr B2 envelope of D/H=(0.2-ll) x 10~® . 
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1. Introduction 

The abundance of deuterium in the interstellar medium is 
a fundamental question in astronomy. Deuterium was only 
produced in significant amounts in the big bang and is 
destroyed via processing in stellar interiors. Observations 
of the D/H ratio in the ISM provide lower limits to the 
primordial deuterium abundance. 

The Galactic Centre is the region of our galaxy with 
the largest star formation rate and contains the most 
highly processed material. Recent mod els of galactic chem- 
ical evolution (e.g. Matteucci et al. 1999 ) predict that 
the abundance of deuterium should be drastically reduced 
in the Galactic Centre due to the large stellar process- 
ing rates. However, recent observations of the deuterated 
molecules DCO + and DCN by Jacq et al. ( |1999| ) and 
Lubowich et al. ( 200C| ) towards Galactic Centre molecular 
clouds show that deuterium is clearly present and not as 
depleted as the models suggest. The final determination 
of the exact amount of deuterium from these observations 
depends both on the line modelling and on the deuter- 



ation process for HCO"'" and HCN. Deuterium in these 
molecules is expected to be enhanced over the actual D/H 
ratio by fractionation effects which increase for lower tem- 
peratures. 

For this reason, a very good way to measure the deu- 
terium abundance is by observing the HD ground state 
rotational transition at 112 /xm. Observation of this line 
measures the HD column density and comparison with 
the relevant H2 column density leads to a straightforward 
determination of the D/H ratio. 

We report on observations of the 112 fim HD line to- 
wards Sagittarius B2 (Sgr B2), a giant molecular cloud 
complex located ~100 pc from the Galactic Centre. The 
complex consists of several clusters of compact Hii re- 
gions and dense molecular cores surrounded by a diffuse 
envelope. Its far-infrared spectrum is dominated by con- 
tinuum emission from dust in the central clusters, with a 
colour temperature of about 35 K and a blackbody peak 
near 90 /im (eg. see the recent review by Ceccarelli et al. 
^0021) . 
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Sgr B2 was observed as part of a wide spectral survey us- 
ing the ISO Long Wavelength Spectrometer (LWS) Fabry- 
Perot (FP) mode (Clegg et al. |l99|). The whole LWS 
spectral range, from 47 to 196 /im, was covered using 36 
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separate observations with a spectral resolution of 30- 
40 km s~^. In each observation both the LWS FP and 
grating were scanned to cover a wide range in wavelength. 
The first description of the survey and its results on am- 
monia lines are reported in Ceccarelli et al. ( ^002 ). 

The LWS beam had an effective diameter of 78" at 



112 ^m (Gry et al. 2002). The beam was centred at co- 
ordinates a = 17'M7'"21.75", 5 = -28°23'14.1" (J2000). 
This gave the beam centre an offset of 21.5" from the main 
far-infrared peak, which occurs near the radio and mm 



source, Sgr B2 (M) (Goldsmith et al. |l99g). This pointing 
was used to exclude the source Sgr B2 (N) from the beam. 

During each observation the LWS FP and grating set- 
tings were optimised for the detector whose band pass 
filter included the wavelength range of interest. This was 
designated as the 'prime' detector. However, all ten LWS 
detectors recorded data simultaneously in their own spec- 
tral ranges. The other nine detectors are known as 'non- 
prime' and often recorded useful data that can comple- 
ment the prime data. 

The region around the HD J = 1 ^ rotational tran- 
siti on at a rest wavelength of 112.0725 /im (Evenson et 
al. 



was covered in two separate observations using 
the LWS long wavelength FP (FPL) and detector LW2. 
The first was a prime observation made on 1997 April 
5 (ISO revolution 506) and the second was a non-prime 
observation made on 1997 April 8 (ISO revolution 509). 
The non-prime observation had limited wavelength cover- 
age, starting at 112.064 /^m (velocity of -23 km s'^). The 
sampling interval used in both observations was a quarter 
of a spectral resolution element with each data point re- 
peated three times. The resolution element, determined by 
interpolating between measurements made on the ground 
and in orbit, was found to be ~33 km s^^. 

Each observation was processed to the Standard 
Processed Data (SPD) level using the LWS pipeline ver- 
sion 8. This produced data still in engineering units (FP 
gap voltage and detector photocurrent). Further calibra- 
tion was carried out interactively using routines from 
the L WS Interactive Analysis (LIA) package (Gry et al. 
2002). The standard LIA tool for processing FP data. 



'FP_PROC', was modified using improved algorithms and 
calibration files in order to calibrate the data in flux (these 
appeared as part of LIA version 10). Due to a large stray- 
light contribution (caused by the strength of the source), 
the dark signal was determined by an algorithm that ex- 
amined instances during on-source observations when the 
grating and FP were aligned to cut out transmission in the 
main beam. The final dark value is ^^3% of the continuum 
photocurrent and so is not a significant source of error in 
the calibrated flux. The relative grating response in FP 
mode was recovered from the non-prime observations in 
the dataset and a correction was applied for contamina- 
tion from adjacent FP spectral orders (for more details of 
these corrections, see Polehampton et al. 2001; Gry et al. 



The two observations were carefully checked for 
glitches and transient affected scans and these were re- 
moved using the ISO Spectral Analysis Package (ISAP). 
The observations were matched in wavelength by correct- 
ing each one to the local standard of rest. The final wave- 
length alignment was determined by examining the posi- 
tion of the H2O 113.540 yum line which occurred in both 
observations. The remaining uncertainty is 0.004 /ini (or 
11 km s~^) which corresponds to the error in absolute 



wavelength calibration (Gry et al. 2002). A small scal- 
ing factor (10%) was applied to the flux in the non-prime 
observation to bring the continuum level into agreement 
with the prime observation. The two observations were 
then co-added. 

The largest errors in the flux at 112 /im were due to the 
extended nature of the source and its offset from the LWS 
optical axis. These can cause a large scale fractu ring and 
bending of the spectrum on each detector (Lloyd 2001 ) . A 
correction to account for these effects was not applied. To 
give an idea of the error on the flux value, the data were 
compared with a medium resolution LWS grating mode 
observation of Sgr B2 taken with the same pointing. After 
correction for saturation in the long wavelength detectors, 
this shows a flux level at 112 /im that is '^30% smaller than 
the FP mode flux level and this reflects the magnitude of 
the error in absolute flux. Other calibration effects in FP 
mode for which the errors can be quantified contribute a 
systematic error of <5%. 

The region around the HD J = 2 ^ 1 rotational tran- 



2002. Further information will appear in Baluteau et al. 



in preparation). 



sition at a rest wavelength of 56.230 /im (Ulivi et al. 1991) 
was also observed as part of the spectral survey. This line 
was covered in three non-prime observations using FPL 
with detector SW2. These observations were processed 
and co-added in a similar way to the 112 /im observa- 
tions. The systematic error in absolute flux is dominated 
by uncertainty in the dark signal which makes up a large 
fraction (40%) of the continuum photocurrent. The dark 
could not be determined in the same way as at 112 /im and 
so an estimate was made by comparison with FP prime ob- 
servations and the medium resolution LWS grating mode 
observation. The final error in fiux is estimated to be not 
more than 30%. 

3. Results 

Figure |l| shows the 112 /im data binned at half a spec- 
tral resolution element. Statistical error bars are shown 
for each bin and these are much smaller than the over- 
all systematic error in the flux (30%). The signal to noise 
clearly increases at velocities above -f 25 km s^"'^ due to the 
combined effect of overlap between FP mini-scans and the 
co-addition of non-prime data. A probable emission line 
detection occurs near the expected velocity of Sgr B2 (M) 
(-1-60-65 km s~^; eg. Mehringer et al. 1992). In order to 
determine if this feature was real, individual mini-scans 
were examined from both observations as edge effects on 
FP scans can sometimes cause spurious features in the fi- 
nal spectrum. However, emission was present on all of the 
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Sgr B2 

HD 1 1 2,0725 ^rr 




Fig. 1. LWS high resohition Fabry-Perot data in the re- 
gion of the HD 112 /im line, plotted against local standard 
of rest (LSR) velocity. The data have been binned at half 
a spectral resolution element. The error in each bin due 
to the spread in the data points is shown. Systematic er- 
rors in flux are not shown and are equal to roughly 30%. 
The best fit to the line is shown after taking account of 
different offsets in the binning grid. The baseline has not 
been corrected for high frequency fringing (see text). The 
expected position of a line from Sgr B2 at -1-60-65 km s^^ 
is indicated. 



overlapping mini-scans before co-addition, showing that 
it is not due to scan edges. The co-added spectrum was 
checked for the high frequency fringi ng th at is sometimes 
seen in LWS FP data (see Gry et al. |2002| ). The standard 
LIA routine, 'Defringe', was used to determine the level of 
fringing, which could not account for the detected feature 
(this routine solves for the period, phase and amplitude 
of a sine wave fitted to the data). The fringing was not 
removed as it only affected the height of the detection 
by 5%. The detection cannot be due leakage from adja- 
cent FP spectral orders (see Gry et al. 2002) as there are 
no strong emission lines present in the nearby continuum. 
The 112 /im line has previously been detected by the LWS 
both in emission towards warm PDR regions in the Orion 



uncertainty in the fit). Systematic uncertainty contributes 
a further error of ~30%. 

The HD J = 2 — > 1 transition at 56 ^m was not de- 
tected in the co-added spectrum. An upper limit for emis- 
sion in this line shows that the line flux must be less than 
1.7x10"^^ erg s~^ cm~^. 

4. D/H Ratio 

Several studies on the dust and gas content of the diffuse 
cloud overlying Sgr B2 (M) have been carried out. The 
structure of the regions around the source have recently 
been summarised by Ceccarelli et al. ( ^002 ). Emission 
from HD could occur in two regions of the Sgr B2 complex 
which appear extended in the LWS beam. These are the 
warm (40-100 K) envelope, with n(H2)« 2 x 10^ cm~^, 
and the hot (~700 K) layer seen in NH3 absorption with 
n(H2)« 10^ cm-3. 

At densities of ~10^ cm""^, the deviation from LTE of 
the HD J = 1 and J = 2 level populations is very small 
(e.g. Bertoldi et al. |1999| ). This means that the total HD 
column density can be derived from the line emission, />,, 
the source extent, risourco, and the gas temperature, Tgas, 



7V(HD) 



47r h 



giexp(-i;i/fcrg, 



(1) 



where Ei/k — 128 K, Ai^ is the Einstein coefficient for 
spontaneous emission, equal to 5.12 x 10^^ s^^ (Abgrall, 



Roucff & Viala p^8^ , gi=3 and Q(Tgas) is the partition 
function at the gas temperature. The angular size of the 
100 /im continuum emission from the extende d env elope 
is approximately 55" x 120" (Goldsmith et al. |l992|) , the 
centre of which was offset from the LWS optical axis by 
21.5". The solid angle covered is « 7 x 10^* sr. 

The observed flux at 112 /im can be used with the de- 
tection limit for the 56 /im line to constrain the tempera- 
ture in the emitting region, assuming LTE. Comparing the 
population in the J = 1 level with the upper limit for the 
population in the J ~ 2 level leads to a temperature in the 
emitting region of Tgas < 80 K (taking account of errors). 
This excludes the possibility that the observed HD emis- 
sion comes from the hot layer (it is unlikely that the hot 
region has a density much lower than 10"* cm~'^, in which 
case the J = 2 level would be subthermally populated). 

Assuming, therefore, that the HD emission occurs in 
the warm envelope (with parameters from above, 7i(H2) « 
2 X 10^ cm-3 and T=40-80 K) gives a HD column den- 
sity ranging from 1.8 x 10^* cm^^ (for r=:80 K) to 
6.0 X 10^® cm~2 (2-^40 K). Taking account of la sta- 



Bar (Wright et al. 1999) and in absorption towards the 
massive molecular cloud W49N (Caux et al. 2002). 

The line width of the absorption component due to 
Sgr B2 in the 6 cm transition of H2C O has been observed 
to be «12 km s~^ (Mehringer et al. 1995 ), which is un- 
resolved by the LWS. Convolution of a Gaussian profile 
of this width and the LWS resolution element at 112 /im 
gives an expected line width Ri37 km s~^. Figure 1 shows 
the best fit of this convolved profile to the data. The ef- 
fect of changing the position of the bin grid over the data 
was investigated and this did not change the fit parame- 
ters significantly. The best fit line is centred at a velocity 
of 51 km s^^ (consistent with Sgr B2 within the LWS (1981). From these observations they derived A^(H2) be- 
FP wavelength calibration accuracy), with a line flux of tween 4 and 5 xlO^'^ cm~^. This estimate is consistent 
(3.3± 1.4) X 10~^^ erg s^^ cm^^ (the error is Icr statistical with studies of the dust extinction towards the central 



tistical and 30% systematic error in the 112 /im line flux 
gives iV(HD) = (0.7-11) x 10^* cm'^. 

Estimates of the molecular gas towards Sgr B2 
have been obtained from studying absorption by trace 
molecules such as H^^CO+ and H^^CN by Link et al. 
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sources (Harvey et al. 1977; Erickson et al. 1977; Gatley 
et al. |1978| ; Thronson & Harper |1986[ Lis et al. |l99lD 



giving A^(H2) = (0.5-2) x 10 cm and observations 
of the whole Sgr B2 complex in sub-mm dust continuum 
(Goldsmith et al |1987| Lis et al. [1993| ; Gordon et al. |1993 ; 
Kuan et al. [1996D . 

Adopting an H2 column density of (0.5-2) x 10^'* cm~^ 
gives a deuterium abundance of (0.2-11) x 10~^, corre- 
sponding to temperatures in the emitting region of 80- 
40 K. These values may still be an underestimate, due to 
the possibility of self-absorption of HD within the Sgr B2 
complex. Evaluating an accurate correction for self ab- 
sorption would require a detailed knowledge of the vari- 
ous components in the Sgr B2 complex and is beyond the 
scope of this paper. However, approximate calculations 
show that absorption is unlikely to change the result by 
more than 50%. 



5. Conclusions 

We report a probable emission line detection of the HD 
J = 1 — > transition at 112 /im towards Sgr B2. An up- 
per limit for the J = 2 ^ 1 transition at 56 /im is used to 
show that the emission must originate in the warm part 
of the Sgr B2 envelope, at a temperature lower than ap- 
proximately 80 K. This leads to a deuterium abundance 
in the Sgr B2 complex of (0.2-11) xlO~^. This range is 
constrained by the temperature (80-40 K) and H2 col- 
umn density ((0.5-2) x 10^** cm^'^) in the emitting region. 
The error in D/H at a particular temperature and A^(H2) 
within this range is approximately a factor of two. 

The D/H ratio derived from DCO+ and DCN obser- 
vations of Sgr B2 by Jacq et al. ( |1999 ) was calculated to 
be in the range (0.35-10) x 10~^, with a firm lower limit of 
10^^. Our values agree with this range and are consistent 
with the deuterium abundance of (1.7±0.3) x 10~^ , foun d 
in the Sgr A 50 km s^^ cloud by Lubowich et al. (2000). 



The chemical evolution o f the Galactic Bulge has been 
studied by Matteucci et al. (1999). They predict that the 
present time D/H ratio in the bulge should be about three 
orders of magnitude smaller than in the Galactic Disk 
due to the shorter timescales for astration processes there. 
Our values show that D/H in Sgr B2 is possibly smaller 
than the average value found in the local interstellar cloud 
((1.5±0.10)xl0"^; Linsky |1998|) but cannot be as low as 



Matteucci et al. ( 1999 ) predict. This may indicate that the 
deuterium content of the bulge has been enriched by infall 
of primordial material or that there could be a problem 
with the model. Prantzos (199(;) has calculated models 
that include infall and shows that the predicted D /H gra- 
dient is very sensitive to the assumed infall history. 

The HD ground state rotational line provides a very re- 
liable way to measure the deuterium abundance as (HD) 
directly traces the deuterium column density. The remain- 
ing errors in the deuterium abundance are mostly due to 
uncertainties in the determination of the H2 column den- 
sity and temperature in the emitting region. Our probable 
detection of the 112 /xm line indicates that with the higher 



sensitivity of SOFIA and maybe Herschel, many more ob- 
servations of HD rotational lines will be possible. An im- 
provement to the deuterium abundance derived from the 
HD line will also require an improvement in the determi- 
nation of the H2 column density. This could be carried 
out either by direct observation of the H2 transitions in 
the near-infrared or by improving modelling of molecular 
abundances. 
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